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rOREWORD 


This  report  was  prepared  by  the  Air  Force  Aero  Propulsion  Laboratory  of  the  Research 
and  Technology  Division  of  the  Air  Force  Systems  Command.  Dr.  Henry  R.  Velkoff  served 
as  project  scientist.  The  work  reported  herein  was  accomplished  under  Project  3141, 
“Electric  Propulsion  Technology.” 

This  work  is  a  continuation  of  a  sequence  of  programs  aimed  at  exploring  tlie  basic 
nature  and  practical  applications  of  the  interactions  of  electrostatic  fields  with  fluids. 

The  general  area  has  been  entitled.  “Electrofluld-mechanics.”  Prior  research  In  this 
area  has  included  (a)  a  basic  review  of  possible  phenomena,  reported  in  ASD-TR-61-642: 
(b)  a  study  of  electrostatic  effects  on  free  convection,  reported  in  ASD-TDR-62-650;  (c) 
a  study  of  electrostatic  effects  on  condensation  of  vapors,  reported  In  RTD-TDR-63-4008; 
and  (d)  an  exploratory  study  of  the  effects  of  ions  on  laminar  flow,  repoj  '^ed  in  ASD-TDR- 
63-842,  This  report  seeks  to  explain  a  phenomenon  observed  in  the  latter  work,  ASD-TDR- 
63-842.  This  research  was  accomplished  from  May  1963  to  September  !''63. 

The  author  wishes  to  express  his  appreciation  to  those  of  the  Aero  Propulsion  Labora¬ 
tory  who  aided  In  this  program,  panicularly  Mr.  Michael  Pennucci  for  his  steadfast 
assistance  In  performing  the  many  computations  necessary  in  this  analysts. 
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ABSTRACT 


An  analysis  was  performad  on  the  actions  of  ions  In  a  laminar  flow  stream.  The  anal¬ 
ysis  seeks  to  provide  an  explanation  for  the  unusually  high  p'essure  drops  found  during 
a  preceding  cxix'rimental  program  involving  laminar  flow  of  ionized  gas  with  an  applied 
electrostatic  field.  The  analysis  revealed  that  very  high  pressure  drops  could  be  ex- 
ivcted  and  that  velocity  profile  distortions  would  occur.  Good  correlation  with  test  data 
was  achieved.  A  new  nondimensional  parameter,  called  the  charge  number,  was  introduced. 
Charge  numlier  is  related  to  the  Knudsen  number  and  is  the  controlling  parameter  for  ion 
action  in  the  gas  stream.  The  conclusion  from  this  work  is  that  ions  bound  electrically  to 
a  flow  stream  set  up  an  internal  field  on  the  stream,  and  then  tend  to  resist  the  motion  of 
the  stream.  This  action  causes  profile  distortions  and  increases  effective  fluid  friction. 


This  report  has  been  reviewed  and  is  approved. 


R'OBERT  E.  SUPP 

Asst.  Chief,  Aerospace  Power  Division 
Air  Force  Aero  Propulsion  Laboratory 
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INTRODUCTION 


During  cn  exploratory  investigation  into  the  effects  of  ionization  on  the  flow  of  gases 
in  a  channel,  unusually  high  pressure  drops  were  observed  when  ionization  was  present 
in  the  channel.  Complete  details  of  that  investigation  may  be  found  in  Reference  1.  In 
that  work,  gases,  such  as  air,  were  passed  through  a  1-1/4-lnch-dlamctcr  pipe  in  which 
a  0.004-lnch-diameter  wire  was  located  concentrically.  A  high  voltage  applied  to  the  wire 
resulted  in  a  corona  discharge  that  provided  ions  in  the  gas  stream.  Under  the  action  of  the 
field,  pressure  drops  in  the  pipe  were  doubled,  velocity  profiles  were  distorted,  and  heat 
transfer  was  doubled.  Typical  data  shown  in  Figures  1.  2,  and  3  Illustrate  the  type  of 
changes  observed  with  the  ion-gas  interaction.  Figure  1  shows  that  data  for  a  given  value 
of  current  fall  on  line.a  that  are  parallel  to,  but  displaced  above,  the  theoretical  54 
line.  This  Indicates  that  the  flow  tends  to  retain  its  laminar  character  even  under  the 
influence  of  the  field.  Figure  1  also  indicates  that  little  influence  is  exerted  in  the  turbu¬ 
lent  flow  regime. 

The  velocity  profile  data  of  Figure  2  are  given  in  terms  of  the  total  pressure  as  meas¬ 
ured  with  a  glass  total  head  probe.  The  velocity  profiles  without  field  apolied  illustrate 
the  distortions  due  to  the  presence  of  the  wire  alone.  vVhen  the  ions  traverse  the  stream, 
very  large  changes  in  the  shape  of  the  profiles  are  evident.  Although  the  data  presented 
are  not  precise  in  the  vicinity  of  the  wall,  the  general  pattern  indicates  a  flattening  of  the 
profile  and  a  corresponding  steepening  of  the  gradient  at  the  wall.  This  action  is  in  the 
correct  direction  to  cause  increases  in  pressure  drop,  or  ,  as  shown  in  Figure  1. 

Data  on  heat  transfer  shown  in  figure  3  provides  further  evidence  of  the  effect  of  ions 
on  fluid  behavior. 

Several  possible  mechanisms  that  could  account  for  the  phenomenon  are  discussed  in 
Reference  1.  It  was  concluded  that  none  of  those  mechanisms  accurately  described  the 
phenomenon,  nor  did  they  fit  the  data  presented.  This  repon  puts  forth  a  hypothesis 
which  attempts  to  explain  the  nature  of  the  phenomenon.  A  series  of  models  based  upon 
selected  approximations  is  presented  to  describe  the  pressure  drop  increase  and  to 
indicate  possible  effects  on  the  velocity  profile.  Comparisons  are  then  drawn  between 
the  models  and  the  data. 


This  repon  was  released  by  the  author  on  5  September  1963  for  publication  as  an  RTD 
Technical  Documentary  Report. 
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Two  approaches  will  be  taken  In  considering  the  actions  of  the  Ions  on  the  stream.  The 
first  is  based  upon  physical  insight  to  the  problem  and  the  second  is  a  more  formal  mathe- 
niathical  approach. 

ION-FLOW  INTERACTION 

In  considering  the  physical  aspects  of  the  problem,  it  Is  helpful  to  utilize  several  highly 
idealized  models  as  well  as  the  concentric  tube  actually  used  in  the  tests.  As  shown  in  I 

Figure  4,  the  gas  flows  down  the  channel  in  the  t  direction.  The  electrical  field  Is  applied  I 

between  the  wire  and  the  tube  wall.  Ions  created  by  the  corona  discharge  drift  from  the 
wire  to  the  wall  under  the  action  of  the  field,  and  at  the  same  time  tend  to  be  dragged 
downstream  with  the  flow.  The  order  of  magnitude  of  the  velocity  normal  to  the  flow  can 
be  obtained  from  the  mobility  equation  for  ions  in  a  gas  (Ref.  2): 

Vj  =  KE,  (1) 

where  V|  is  the  ion  velocity  normal  to  the  stream, 

K  is  the  mobility  of  the  ion.  and 

E,  is  the  electrical  field  strength. 

For  the  case  under  consideration,  typical  values  are: 

2. 0  X  to'*  w* 
voii-t#c 

and 

E  s  5  X  10*  voM»/m. 

Thus,  Vj  is  of  the  order  of  100  m/sec. 

In  the  case  of  laminar  flow  with  air  at  standard  conditions,  wm  ■  1  m/sec,  where 
is  the  mean  stream  velocity.  If  the  Ion  is  dragged  downstream  with  the  flow,  it  will 
acquire  approximately  this  value  of  axial  drift  velocity.  The  net  effect  of  this  condition 
is  that  the  ion  travels  almost  directly  across  the  channel  under  the  influence  of  the  applied 
field. 

If  the  above  description  is  assumed  valid,  then  it  does  not  appear  possible  for  the  ions 
to  interact  with  the  flow  stream.  Another  interpretation,  however,  is  also  possible.  If  the 
Ions  were  to  drift  at  all,  then  this  would  imply  that  some  of  tliem  would  be  carried  oat  of 
the  pipe  at  the  end.  Such  an  action  would  violate  the  condition  of  current  continuity,  that 
is,  the  ions  leaving  the  wire  must  collect  on  the  outer  electrode  in  order  to  complete  the 
electrical  circuit.  If  this  condition  is  met.  then  the  ions  must  try  to  stay  within  the  confines 
of  the  pipe.  To  do  so,  some  type  of  Internal  field  must  exist  to  cause  the  Ions  to  move  in 
opposition  to  the  flow.  As  a  first  order  approximation,  it  can  be  assumed  that  the  ions 
resist  being  dragged  by  the  axial  velocity  of  the  stream,  and  as  a  limiting  case  would  tend  | 

to  move  directly  across  the  stream  without  any  net  w  velocity.  For  this  limiting  case,  the 
ions,  on  the  average,  would  tend  to  remain  stationary  relative  to  a  given  position  along  the 
pipe  as  they  move  In  the  radial  direction  under  the  action  of  the  applied  field.  The  velocity  j 

of  the  ion  is  given  by: 
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where  J  is  the  stream  velocity  In  jl!!.  direction, 

uj  Is  the  Ion  velocity  In  the  j  direction,  and 
u^ls  the  ion  velocity  relative  to  the  stream  velocity. 

The  mobility  equation  for  an  Ion,  In  general  form,  is: 

4  *  ‘‘e' 

and  in  the  radial  direction  is: 

V;  =  V  4- 

=  V  +  KE,  . 

But  the  average  stream  velocity  In  the  r  direction,  v ,  is  zero.  Hence, 

V|  KEf  , 

which  Is  Equation  (I).  In  the  case  of  the  axial  motion 

Wj  e  w  4- 

wj  *  w  4-  KEj  . 


(2) 


(3) 


(4) 


(5) 

(6) 


If  the  ion  drifts  with  the  stream,  then  is  zero  and  w  :  wf  ,  and  no  field  exists.  If, 
however,  the  other  limiting  case  applies,  w|  goes  to  zero  since  the  ions  remain  stationary 
relative  to  the  walls.  In  this  case: 


•i 


—  0 


w  *  —  KE* 

This  equation  indicates  that  an  internal  field  may  be  established  in  the  gas  stream  which 
is  a  function  of  the  stream  velocity  and  mobility: 


As  a  limiting  case. 


E* 


(7) 

(8) 


An  estimate  of  the  order  of  magnitude  of  the  electrical  fields  Involved  can  be  determined 
by  using  Equation  (8)  and  assuming  typical  values  for  K  and  w, 

-4  i  / 

K  *2. ox  10  m  /volf-8*c  ond  w*  lm/s«c. 

Then: 

“  T 

E  I  *  —  50  volf  t  /cm. 
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Thus,  the  Induced  fields  are  small  relative  to  the  applied  radial  fields. 

If  the  usual  assumptions  for  laminar  Hagen-Poiseuelle  flow  are  made,  the  Navler-S'.okes 
equations  become: 


^  V  *  «•  s  — W 
di 

where  p  is  the  pressure,  is  the  body  force  In  the  z  direction,  and/i  is  the  viscosity  of 
the  gas.  Assuming  the  only  active  body  force  is  due  to  charge  action  in  the  stream, 

^Z  ^Pc  ^z 

where  is  is  the  scalar  charge  density: 

+Pc  Ez  *  0-  (") 

Substituting  Equation  (8), 

-p  ^  .  (12) 

®  dz 

This  equation  can  be  solved  once  a  charge  density  distribution  is  known  using  the  usual 
assumption  that  dp/dz  is  a  constant  for  fully  established  steady-state  pipe  flow.  This 
equation  represents  the  limiting  case  for  w;  z  o  . 

Before  considering  solutions  to  this  equation,  a  more  extensive  consideration  will  be 
given  to  the  effect  of  current  continuity  on  the  ion  motion  phenomenon.  The  equations 
applicable  to  the  problem  at  hand  are  a  combination  of  electric-field  and  fluid-flow  equa¬ 
tions  (Ref,  3).  The  problem  attacked  here  is  that  of  the  isothermal  flow  case;  the  energy 
equation  is  not  Included  since  the  energy  of  corona  discharge  is  so  small  (Ref.  1). 

ELECTRIC-FIELD  AND  FLUID-FLOW  EQUATIONS 


P  »  P  g  RgT 

(13) 

dp  A 

-  +  — ^  (  PU:  1*0 

dt  d*j  ^  ' 

(14) 

Ou  do  M  A  ^ 

P  * - 4-iiV  0  1  +  +  Fi 

O'  ^  '  3  dKj  ‘ 

(15) 

«  ±L 

(16) 

(17) 

y  X  H  »  7  +  ^ 

dt 

(18) 

„  -  d  u.'h 

V  X  E  = - — - 

dt 

(19) 
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aji 

-  +  -  S  0 

d  \  d  xj 

jj  =  i*  u|  =  a  [ej  +^,  (  u  XH)]  +  u{ 


V*  H  =  0 


V 


E 


u 


j 

R 


=  KE^ 


(20) 

(21) 

(22) 

(23) 

(3) 


ASSUMPTIONS  FOR  CHANNEL  FLOW 

The  following  assumtions  are  made: 

1.  The  stream  velocities  are  very  low,  of  the  order  of  0  to  2  meters/second,  and  M,— O. 

2.  The  fluid  is  Incompressible. 

3.  Flow  Is  steady-state;  no  local  time  derivatives  exist. 

4.  Flow  Is  fully  established  In  the  channel;  the  entrance  region  will  not  be  considered. 

5.  There  are  no  gravitational  body  forces;  In  the  Isothermal  case, this  Is  considered 
a  valid  assumption. 

6.  No  magnetic  effects  are  present.  No  external  magnetic  fields  are  applied,  and  the 
currents  involved  are  so  small  as  to  make  induced  magnetic  fields  negligible.  No 
time-changing  electrical  or  magnetic  fields  are  present. 

7.  The  nonuniform  gradient  terms  are  small  relative  to  any  charge  motion  effects. 

This  is  particularly  true  for  a  nonpolar  gas.  but  It  may  not  be  always  true  for  polar 
liquids  when  highly  nonuniform  gradients  !i»'e  present. 

8.  The  electrical  conductivity,  cr .  of  the  gas  is  extremely  small.  It  is  assumed  that 
all  current  is  due  to  the  ion  motion. 

With  these  assumptions,  the  field  and  flow  equations  become: 


3ui 


3*1 


/>  u  — 


(24) 

(”5) 
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'^j  '  Pc 

(20) 

7  X  E  s  0 

(27) 

dt  dx^ 

(28) 

=  o-E^  +  /),  u| 

(29) 

7-  E  =  -5e- 
€ 

(30) 

(3) 

Considering  liquations  (30)  and  (26): 

*  *  7  •  F 
'  c 

Fj  »  Ej  *  €  Ej  7  •  F 

'i  •  (-ffr)  ■  P‘) 

For  cartesian  coordinates: 


6 

dE, 

dEy 

dEy  »1 

s  1 

dx 

•t* 

dy 

' 

di  /  ^ 

A 

+  i 

[«Ey  1 

f  dE, 

dEy 

dEy  V  1 

^  dx 

▼ 

dy 

"r 

dz  '  J 

A 
+  k 

[‘E,  ( 

’  dE, 

dEy 

X 

'  d  X 

T 

d  y 

T 

dz  '  J 

Thus,  the  body  force  can  be  expressed  explicitly  if  cither  the  charge  density  distribution 
or  the  field  is  known.  In  rlie  case  of  two-dimensional  flow  between  flat  plates  with  a  two- 
dimensional  field.  Equation  (32)  becomes: 


jl. 

dZi 

-4-  *  Ey 

dEy 

(33) 

dy 

dz 

♦E, 

dE, 

+  -i-  - 

UJ 

(34) 

dy 

2 

dz 

The  terms  on  the  trace  represent  the  gradient  of  the  ordinary  electrostatic  pressure  rise. 
The  other  terms  indicate  the  possible  cross  influences.  In  particular,  if  any  Ej  field  at 
all  exists  in  the  flow,  the  dz^/d^  term  may  be  large  enough  to  be  significant. 
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With  no  sources  of  Ions  throughout  ihe  stream,  a  quasl-steady-state  process  exists. 
Then  Equation  (28)  becomes; 

»  O  .  (28a 

d.' 

For  example,  consider  the  parallel  plate  case  again.  Using  Equations  (28a)  and  (29)  with 

o-  =0, 

dp.  u! 

- -  *  7-  *  *  0 .  (35) 

dx:  OU: 


For  the  two-dimensional  case: 

.  o .  (36) 

3  y  di 

As  an  analogous  case  to  the  flow  of  Ions  from  tlie  wire  to  the  wall,  assume  that  the  corona 
discharge  takes  place  at  the  flat  channel  wall  and  the  ions  drift  across  the  space  to  the 
other  wall.  In  this  case,  the  net  Ion  current,  | ,  will  be  given  by: 

1  =  J^A  (37) 

where  a  the  cross  section  of  the  flat  plate  over  which  the  fluid  flows, 

i  =  p^  V.  A  .  (38) 

Although  and  ».  vary  with  »  .  i  Is  a  constant.  Substituting  Equation  (38): 

±  j£sJ^l±  ^  jLPc^  ,  0  (39) 

*  3y  3  / 


Therefore: 


For  fully  established  flow. 


For  p  constant. 


“i 


A.  *i  *  f  <»l 


Pc  "  +  A»c  "r  *  f  <  X » 


Pg  (  w  +  KE  ,  »  =  f  (  y)  . 


w  *  f,  (  y  )  only . 


E*  »  ff  ly )  . 
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Let: 


«»  =  X  On  v" 
nso,i 

(46,1 

E,  =  I  yO  . 

nso,! 

(47) 

From  Equation  (29), 

j'  =  )+crE'. 

(29) 

For  a  — 0 

j*  =  p^  (u^  -t-  )  . 

(48) 

Substituting  Equation  (3): 

j*  --P^{J  +  ke'  )  . 

(49) 

In  they  direction: 

Jy  *  Pf.  ^  Pc  ■ 

(50) 

Since  w  *0  in  channel  flow: 

Jy  */>c  ''b  *  Pc  • 

(51) 

This  is  the  equation  of  current  density  for  the  case  of  corona  discharge  in  a  gas  with 
parallel  electrodes.  In  the  t  direction: 

w  KEj.  (52) 

i'hc  physical  conditions  on  the  test  require  that  no  net  current  flows  out  the  ends  of  pipe: 


I  *  /  dA  »  0 
A”' 

l'lx:rcfore: 

J(^^E^)dA*-  wdA  .  (55) 

Lonsitlcr  two-dimensional  channel  flow  with  a  channel  of  height  2h  ;  Equation  (55)  tiien 
l'-•el>m(.•s: 
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Letting 

y 

V  =  — 

I  I 

h  J  p  C.  ^  d  Tj  -  -  — ^  J  p^  VI  drj  , 

-i  ^  *<  -I  ® 

Rewriting,  Lquations  (46)  anJ  (47)  become: 

w  =  Z  A„  7,"  ;  E,  *  I  B„  7," 

and  letting  lie  constant: 

f.1 

J  1  (A„  4-  K0^  )  7)"  (It;  =  0 

-  I 

I  «A„  +KB„)  -^1  =  O  , 
n*o,i,2 

if  this  function  were  evaluated  over  an  indefinite  interval  instead  of  over  the  specific 
Interval  -1  to  1,  then  such  an  expression  would  lead  to  the  conclusion  that  each  of  the  co¬ 
efficients  would  have  to  be  zero.  If  this  were  the  case,  then: 


Z  A„  7)"  s-I  K  B„  7," 


and 


E,  *- 


Because  of  the  definite  integral,  the  sum  becomes 

An  +K  B 


f.  "  "  )  V  o  . 

'  n  4-  I  f 


With  the  conditions  available.  It  Is  not  possible  to  evaluate  the  coefficients  Ap  In  terms  of 
Bp  .  It  does  not  appear  oosslble  to  uniquely  determine  the  coefficients  of  the  power  series. 

Considering  Equation  (55),  we  can  see  that  Equation  (8), 

satisfies  the  equation  identically  and  consequently  it  is  a  permissible  relationship.  It  Is 
not  likely,  however,  that  It  is  a  unique  relationship. 

The  Navier-Stokes  equations  for  fully  established  flow  become; 

dp 


flV  w  — 


d  X 


AL  Ex  '  0 . 


(11) 


Solving  for  E,  and  substituting  into  Equation  (54),  we  find: 


9 


RTD-TDR-63-4009 


/  - L_  .0, 

A^'  H-  d  i  ' 

One  condition  for  the  flow  which  satisfies  the  current  condition  is: 

I  3p 


V*  «, - ^ 

/tK 


d  2 


s  0  . 


(56) 


(12) 


This  equation  may  he  considered  to  offer  a  first  approximation  to  the  influence  of  the  ions 
on  the  flow  stream.  As  mentioned  above  for  Equation  (8),  however,  this  is  not  necessarily 
a  unique  condition. 

Equation  in5)  can  lead  to  an  insight  into  the  phenomenon  involved.  From  the  definition 
of  the  mean  value  of  a  function: 


-r  ./ 


(5/) 


-i-  r  /)  -ii-dA  *  /-^c  *  .\ 

The  mean  value  of  the  body  force.  Equation  (26),  becomes, 

~  ^  Pc'" • 

*m  K  rj  m 

Substituting  into  Equation  (11), 

dp  I 


«0. 


(58) 


(59) 


(60) 


In  summary,  it  appears  that  ns  a  first  approximation  the  equation  of  motion  for  tlx;  flow 
with  ions  present  in  a  field  can  be  given  by: 


rt  t  Pq  dp 

*  ty  .1  _ 


(12) 


In  addition,  based  upon  the  mean  value  of  the  body  force,  the  equation  of  motion  can  be 
given  by: 


/I  V  *  w - -  (  A.  w  )_  —  — ^  =  O 

K  ^  dt 


(60) 


where  i  ^  »  1^  is  a  constant  value.  Finally,  It  Is  logical  to  assume  that: 
and  Equations  (8)  and  (59)  indicate  two  permissible  forms  of  the  relationship. 
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N 


ANALYSIS 


Equations  (12)  and  (60)  for  the  case  of  channel  flow  In  a  constant  area  duct  with  incom¬ 
pressible  flow  become  ordinary  linear  differential  equations.  Equation  (12)  contains  a 
variable  coefficient,  ,  while  Equation  (60)  has  only  constant  coefficients. 

NONDIMENSIONAL  EQUATIONS 

Before  proceeding  with  the  solutions  to  these  equations,  they  will  be  put  Into  non- 
dimensional  form. 

Let; 


w  = 

w 

• 

(61) 

V  ■ 

r 

R 

for 

tht 

cylindricol  aoso,  and 

(62) 

If 

y 

h 

for 

fht 

porolltl  plot#, 

(63) 

where  w„,  is  the  mean  stream  velocity, 

R  Is  the  Inside  radius  of  the  tube,  and 
h  is  the  half  height  of  the  channel. 

For  channel  flow: 

■ 

and 

V*w  ■ 

FLAT-PLATE  EQUATIONS 

Equation  (12)  becomes: 

d*  W 
d 

The  term  (-4®-)  can  h 

P  •m  dt  2  dC 


a**. 


j —  for  paroliol  plotot 


d  w  .  I  dw  . 

— — *“  ♦ - for  0  round  fubo. 

dr  r  dr 


P 


W  e 


MX 

e  rewritten  as: 


(64) 


(65) 
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n*- 

p 

(66) 

P  • 

P  •'m* 

2 

II 

z 

h 

(67) 

Nr,  = 

h 

(68) 

II 

X 

o 

4A  .w  . . 

-  af  4h 

W.  P. 

(69) 

d*  W 

1 

E 

1 

z 

1:5 

• 

dp* 

(70) 

dT7* 

/J.K  2 

n  ■ 

All  terms  of  Equation  (70)  are  dimensionless. 

Since  w  is  dimensionless,  then  must  be  dimensionless  as  well.  It  represents  the 

unique  contrihutlon  of  the  Ions  on  the  flow,  and  Is  the  characteristic  number  for  this  type 
of  flow. 

Let : 


Therefore 


Define: 


Pc 

Pen.  ■ 


fiK 


(71) 


(72) 


(73) 


where  p  is  the  mean  value  of  charge  and 

L  is  the  characteristic  dimension  of  the  flow. 


This  dimensionless  number  represents  the  ratio  of  the  charge  forces  to  the  viscous  forces 
In  the  flow. 


For  flat  plates.  Equation  (73)  becomes; 


Substituting  into  Equation  (64): 


d*  W 


7-i - N.  p*vt  .  JeI 


(74) 
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This  Indicates  that  tlv:  flow  depends  upon  two  independent  parameters,  the  flow  charge 
number  and  the  Reynolds  numter.  It  should  be  kept  in  mind  that  can  be  a  variable. 

CYLINDRICAL  GEOMETRY  EQUATIONS 


Substituting  into  Equation  (12),  we  find: 


or: 


d*  w 

dT,* 


^  -(-A!L)w 

di7  '  /iK  ' 


_£p_ 

M  "m 


(75) 


d"w 

Tv 


I  dW 

7J  dTJ 


(76) 


where 


P 

M 

P 


c 


*  ■  H- 

'c  ^  K 


(66) 

(77) 

(78) 


The  same  two  independent  dimensionless  parameters  are  Involved  as  before.  The  terms 
in  Equations  (65)  and  (75), 


h 

^  '»m 

are  constant  during  Integration. 


J2_ 

dt 


ond 


EQUATION  WITH  MEAN  VALUE  OF  BODY  FORCE 
Equation  (60)  can  be  put  into  nondimensional  form. 

Flat  Plates 

^  _iEl  (79) 

'  ’  '’'m  ■”  '■  ' 

Round  Tubes 

jIjl  ^  j_  _!«l  N  .Nr  (80) 
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SOLUTION  OF  THE  MEAN  VALUE  EQUATION 

Considerable  insight  into  the  action  of  the  charge  on  the  stream  can  be  gleaned  from  an 
analysis  of  the  approximate  mean  value  equations.  Consider  Equation  (79).  It  can  be  written 
in  the  form: 

d^w  .  h*  /  dp  V  ^ 
dij  ^  di  fiK 

z  ***  /  ^  ”  )m  \ 

\  dt  K  /  ' 


(^c  -L 


Define: 


Therefore: 


Since  the  term  on  the  right-hand  side  Is  a  constant,  the  solution  is  that  of  Hagen- Foiscucllc 
flow: 

W  *  -  — - - n  -17*  )  (82) 

.  :  -  4^n-77*)  (83) 

Zfi.  Az  ' 

•  ,  --Hi 

max  Az 


,  ,  __L  JlL  (1 

">  3  M  Az  • 

Ap 

Substituting  Equation  (81): 

(  h*  /  dp  .  (^c  *)m  \ 

*  T  “  —ir~) 

_  .  _i_  J<1_  (^*)m  ,  _  J _ !L.  ( 

3  M  K  ^  H-  di  ' 

This  equation  indicates  that  for  a  given  mean  flow  velocity,  the  pressure  drop  will  in¬ 
crease  with  the  amount  of  charge  present.  To  more  clearly  show  the  effects,  consider 

the  case  of  constant  charge  density.  Assume  Is  constant.  Then: 


{Pc  *)m  ‘  -m  . 


Equation  (85)  becomes: 
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_L  JiL/_££_\ 

3  M  '  <32  '  /oA^ 

"m  * - 

V  i  fiK  / 

Equation  (S5)  ln(iicatefi  that  for  a  given  available  pressure  drop,  the  flow  velocity  will 
decrease  as  the  charge  density  increases.  This  action  Is  In  agreement  with  the  test  data 
obtained. 

Consider,  next.  Equation  (80)  for  the  round  tube: 


a*w  ^  _l_  dW  _  R*  dp  ^ 

dTj*  7j  di7  di  M  •<  “m 


Defining: 


I  {Pr  ")r 
Lz  di  K 


d'W  .  I 


d  TJ*  ‘T)  a  7} 
The  solution  of  this  equation  leads  to: 


H-  wm  Lt 


-TUL 

Az 

R* 

AP, 

Az 

AP, 

(I-  7)*) 


8  /i  A  I 

Substitution  of  Equation  (87)  leads  to: 

R*  i^c  *)m  ,  _  R^  dp  ( 

8/1  ^  Bfi  dz 

which  is  similar  to  Equation  (85). 

For  the  case  of  constant  charge  density: 

R^  dp 
Bfi  dz 

•m  * - i -  * 

TTaZT 

'  8/1  K  ' 

which  is  similar  to  the  form  of  Equation  (85)  for  flow  between  infinite  parallel  plates. 
Expressing  these  relationships  In  terms  of  the  charge  number,  Np  : 


For  Flat  Plates 


2 


“m 


h 

Ifi. 


dp 

d  Z 


1-^) 


For  Round  Tubes 
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(94) 


R  dp 
Bp.  dz 

The  Blaslus  friction  factor  Is  defined  as  follows: 

T„ 


f 


B 


(95) 


(96) 


where  is  the  fluid  shear  stress  at  the  wall.  The  friction  factor  can  also  be  defined  in 
terms  of  the  pressure  drop  in  the  channel. 


For  Round  Tubes 


»8  ‘2 


«R, 


For  Infinite  Parallel  Flat  Plates 


N 


(97) 


B 


h*  dp 
L.  dt 


Nr,  » 


"m  Nr, 
P  "m  Dh 


(98) 


The  friction  factor  ordinarily  used  in  engineering  work  is  the  Darcy  friction  factor, 
which  is  defined  as: 


'o  •  ^  ^B  •  (99) 

Both  versions  of  friction  factors  arc  used  in  this  report. 

Based  upon  Equations  (94)  and  (95),  the  friction  factors  for  flow  with  charge,  assuming 
p  is  a  constant,  become: 
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CHARGE  DENSITY  IN  THE  CHANNEL 


In  order  to  determine  the  significance  of  the  charge  number  on  the  flow,  It  Is  necessary 
to  determine  t’ne  charge  density.  Its  magnitude,  and  Its  distribution  In  the  channel.  There  are 
two  general  approaches  for  determining  the  distribution  of  charge  of  a  corona  discharge  In 
a  dense  gas. 

LOED  SOLUTION 

Locb  and  Brown  Indicate  a  method  suitable  for  determining  the  distribution  of  charge 
between  coaxial  electrodes  (Refs.  4,  5).  Their  method  assumes  that  with  small  corona 
currents,  the  space  charge  will  be  very  low,  and  the  distortion  of  the  field  lines  will  be 
negligible.  Parts  of  Loeb’s  analysis,  for  coaxial  electrodes,  will  be  reproduced  here. 

The  motion  of  the  ions  Involved  in  the  corona  current  is  assumed  to  act  essentially 
independent  of  the  flow  velocity.  The  problem,  then,  is  reduced  to  determining  the  current 
motion  in  a  radial  direction  only: 


‘  *  2  r  ^ 

where  £  is  the  length  along  the  pipe.  Assume  ■£•  unity.  Substituting  Equation  (3): 


i  s  2  JT  r  p 
*  c 


KE, 


Er  * 


v-T 


-  Zrrp^r 


2  ir  K  V 

in  ^ 


2ir*K 


4  4»r'r*K]  . 


(107) 

(108) 

(109) 

(110) 


For  small  charge  density,  Loeb  assumes  that  the  second  term  can  be  dropped.  Then: 


2ir  K  V 

i  a  — : — ^ 


(111) 


p  a  - ^ -  J^n  -p— 

2f  K  V  "w 


(112) 


Loeb  points  out  that  this  approximation  holds  in  a  region  outside  of  the  intense  corona 
surrounding  the  wire.  The  radius  of  the  intense  corona  region  at  a  pressure  of  one  atmos¬ 
phere  is  of  the  order  of  10^  mean  free  paths.  At  standard  conditions;  this  distance  is  of 
the  order: 


A  r  »  0  <  10*  *  1  cm 


(113) 


a  0.009  inch  >  0.0076  cm 


where  fe  Is  the  radius  at  the  edge  of  the  corona.  In  terms  of  and  a  tube  radius  of  0.625 
inch. 
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7j^  a  0.  005 


Thus,  under  these  assumptions,  the  space  charge  density  will  be  essentially  constant 
between  the  electrodes. 

STUETZER  SOLUTIONS 

Another  approach  for  obtaining  chaise  density  distributions  is  presented  in  summary 
form  by  Stuetzer  in  Reference  6.  The  space  charge  effects  are  included  in  this  approach, 
and  it  is  assumed  that  the  field  at  the  intense  corona  goes  to  zero  because  the  discharge 
is  space-charge  limited.  Only  hip  results  are  presented  here  for  ready  reference. 

Coaxial  Electrodes 


«  .  f  “  .  '  1* 

'c  1  ZrRKi  (r*-rc*)  ^ 

(115) 

"  1  ZirRjtK 

(116) 

Parallel  Plate  Electrodes 

,  .  I  “  f 

>•  2KA(y-y^)  J 

(117) 

.  r  2Hy-yo)  jT 

A  «  K 

(118) 

where  A  is  the  cross  section  normal  to  ion  current  flow,  and  y^ 

is  the  distance  to  the 

edge  of  the  Intense  corona. 

In  the  region  away  from  the  intense  corona,  r  >  and  the  following  approximation 
for  Equation  (115)  results  Jn; 


.  I  (- 

I  RK  J 

where  is  the  current  density  at  the  wall. 


(^r  (-?-)" 


The  average  charge  density  for  a  region  is: 

<  Im  I 


V..  ■/  (t7^)-;T5 


J.  (JH-V 

3  \  aie  I 


t  \«  r'/«  I  */- j 


19 


RTD-TDR-63-4009 


CHARGE  /VND  CURRENT  DENSITY 


\  RK  / 


(122) 


To  calculate  the  actual  magnitude  of  the  charge  densities  used  In  the  tests,  the  actual 
test  configurations  must  be  considered,  and  the  mobilities  of  the  ions  must  be  known. 

Exploratory  Sequence 

i 

R 

Revised  Test  Sequence 

J 

R 

Table  1.  Mobility  of  Ions 
(From  Reference  2) 


I  72  inchts  ^  1.826  maters 
=  0.625  In  =0.01587  malar  t 


a  67  inchas  =1.70  malar  s 
=  0.625  In.  =0.01587  malars 


It  must  be  recognized  that  ion  mobilities  are  not  precise  quantities.  The  values  shown 
are  greatly  affected  by  age,  contaminants,  moisture,  and  to  some  extent  by  the  field 
strength.  As  in  most  corona  wind  analyses,  however,  we  have  assumed  that  K  is  constant 
(Refs.  7,  8). 

For  both  sequences  of  tests  reponed  in  Reference  1,  the  current  density  and  the  charge 
densities  based  upon  Equations  (112)  and  (122)  are  given  by: 
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Exploratory  Teste 


Jw  -5.49  1 

Constant  charge  density.  Equation  (112) 


For  Equation  (122) 


Final  Tests 


Equation  (111): 


Equation  (122): 


fc  -  »■  W 


V  •  ’  *  ^ (i-r 


*  5.88 


Pc  '  ° 


s  7.65  X  10'*  {-^)* 


where  Is  given  In  amp/m  , 

I  is  given  in  amperes, 

3 

is  given  in  coul/m  ,  and 
V  is  given  in  volts. 

To  calculate  the  charge  number  density,  it  is  necessary  to  multiply  the  equations  of 
charge  density  by  the  number  of  charges  per  coulomb: 

I  coulomb  »  1.602  X  lo'*  chorgt*.  (129 

For  a  typical  case,  assume  the  following  values,  based  upon  test  data: 

i  E  I  milliomp 
V  z  6500  volts 
K  *  2  X  lO'^  m*  /  volt  -stc 


For  the  exploratory  tests,  these  conditions  lead  to: 
Equation  (124) 

p^  *  5,92  X  lO*  cborgts 
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Equation  (125) 

p  =  2.56  y  10*  chorgat  /  cc. 

For  this  sat  of  conditions,  values  for  Equations  (115)  and  (119)  are  plotted  in  Figure  6 
along  with  the  average  values  calculated.  It  can  be  seen  that,  except  for  the  region  close 
to  the  wire.  Equation  (119)  gives  a  good  approximation  to  the  space-charge-llmlted  varia¬ 
tion.  Although  the  two  average  values  of  charge  density  are  not  the  same,  they  are  in 
good  agreement  from  an  order-of-magnitude  standpoint. 


One  of  the  parameters  used  in  the  analysis  is: 


(130) 


The  other  factors  used  In  the  analysis  are  related  to  Y^  follows: 
From  Equation  (122) 


From  Equation  (119) 


(131) 


(132) 


From  Equations  (131)  and  (73)  for  the  round  tube  and  variable  charge  density,  and  since 


Pc 


0 


Y, 


MAGNITUDE  OF  THE  INFLUENCE  OF  CHARGE  NUMBER 


(133) 


A  typical  value  for  the  charge  number  can  be  calculated  as  follows.  For  air  at  standard 
conditions: 

K  *2.0  XIO  *  m*  /*oU-s*c 
/i  *  I.  783  X  10'*  kg  /  m  -  $*c 

R  *  1.587  X  10'*  m 


R* 


*  7.  06  X  10* 


m*  / c out. 


From  Equation  (125)  for  the  exploratory  tests  and  a  current  of 


amp: 
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=  7.4  X  io‘*  X  v'T" 

=  16.  5  X  lO’*  coul  /  m* 


Therefore: 


That  this  Is  a  surprisingly  large 
for  friction  factor: 


number  can  be 


II.  7 

seen  from  a  consideration  of  the  equation 


Nr,  ^ 


I  + 


(101a) 


64 

Nr, 


I  2.46) 


(134) 


Assuming  the  charge  density  given  by  Equation  (124)  and  I  >  lO**  amp  leads  to: 


(135) 


Thus,  the  magnitude  of  the  charge  number  is  sufficient  to  cause  increases  in  friction 
factor  which  are  similar  to  those  found  in  the  test  data.  From  consideration  of  Equations 
(73),  (101),  and  (125),  It  can  be  seen  that  the  Increase  in  friction  factor  depends  on  the 
square  root  of  the  corona  current: 


64 

Nr, 


I 

I  1+  oi^) 


(136) 


where  a  Is  a  constant. 

Equations  (100),  (101),  and  (134)  predict  that  tie  curves  for  friction  factor  should  be 

parallel  to  the  64/ Nr  line  for  flow  without  Ionization.  Values  computed  from  these 

• 

equations  and  the  test  results  plotted  In  Figure  1  are  strikingly  similar.  Aithough  the 
absolute  values  of  fQ  as  given  by  Equation  (134)  do  not  agree  precisely  with  the  test  data, 
the  degree  of  correlation  achieved  is  surprlsin^y  good.  It  must  be  borne  in  mind  that 
these  expressions  for  friction  factor  are  based  upon  a  mean  value  approximation  to  the 
body  force  which  was  used  in  the  solution  of  the  differential  equation  and  that  the  mean 
value  approximation  leaves  the  shape  of  the  velocity  profile  unchanged  as  the  ions  flow 
throug*'  the  gcs.  Test  data,  however,  show  a  marked  influence  on  the  velocity  profiles. 
Consequently,  to  evaluate  the  possible  effects  of  the  charge  on  the  profiles  as  well  as  on 
the  pressure  drops,  it  is  necessary  to  consider  the  equations  in  which  the  body  force 
varies. 
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SOLUTIONS  FOR  CHANNEL  FLOW  WITH  CHARGE  PRESENT 


Six  separate  cases  for  flow  through  a  channel  are  considered.  They  include  variations 
in  the  flow  boundary  conditions  and  two  different  assumptions  of  charge  density  distribu¬ 
tion,  as  follows: 

Case  1.  Infinite  parallel  plates  with  constant  charge  density,  p.  . 

*^0 

Case  II.  Infinite  parallel  plates  with  variable  charge  density  given  by  Equation  (117). 

Case  111.  Pound  tube  geometry  with  constant  charge  density,  ,  assuming  the  wire 
does  not  influence  the  flow.  ° 

Case  IV.  Round  tube  geometry  with  constant  charge  density,  p^  ,  assuming  the  flow 
Is  at  rest  at  or  near  the  wire.  “ 

Case  V.  Round  tube  geometry  with  variable  space  charge,  assuming  the  wire  does  not 
affect  the  flow. 

Case  VI.  Round  tube  geometry  with  variable  space  charge,  assuming  the  flow  is  at 
rest  at  or  near  the  wire. 

CASE  I.  PARALI  EL  PLATES,  A  CONSTANT 

c 

This  is  the  simplest  idealized  case  to  consider.  We  have  assumed  that  a  uniform  charge 
exists  between  two  infinite  paraUel  electrodes.  Equations  (64)  or  (74)  are  applicable. 

Consider  Equation  (74); 


d*  w 


nr. 


d  7) 


-  -N«  r 

t  fc  2 


Since 


for  this  case: 


»♦  - 

«  Pc 

'•m 


»  I 


d  W  Nr, 

d,,*  ■  ‘  2 


ac 


_h 
M-m 


d  *  • 


This  equation  can  be  solved  directly  with  the  boundary  conditions  of. 


y  * 

0, 

dW 

4y 

‘0, 

• 

y* 

h. 

w  = 

0  . 

w* 

_  / 

f  h*  ' 

\  dp 

1 

'  di 

CO»h  ./Npe  Tf  1 

cosh  >/N  p^ 


(74) 


(137) 


(138) 

(139) 

(140) 
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As  «  —  0  .  this  expression  becomes: 

* C 


“-T  — 

2  uw-  \  di  f 


which  is  the  equation  for  flow  without  ionization. 

It  Is  of  interest  to  compare  the  form  of  Equation  (140)  with  the  Hartmann  equation  of 
magnetohydrodynamics.  From  Reference  9: 


f  ^  1  Ui 

coth  M  1 

cosh  M  * 

where  Is  the  a.opHcd  magnetic  field. 

Is  the  permeability  of  the  medium,  and 
M  is  the  Hartmann  number. 


The  term  in  the  left  set  of  brackets  of  Equation  (141)  is  a  constant  for  a  given  set  of 
parameters.  Likewise,  the  terms  preceeding  tlie  bracket  of  Equation  (140)  are  constant. 
We  can  readily  see  that  Che  two  profiles  are  of  the  same  general  form.  In  MHD,  this  flow 
leads  to  a  flattening  of  the  velocity  profiles.  It  can  be  expected,  therefore,  that  Equation 
(140),  which  deals  with  charge  motion,  leads  to  a  similar  flattening.  This  flattetilng  is 
apparent  in  Figure  7,  where  two  values  of  the  charge  parameter,  N*  ,  for  Equation  (140) 

"c 

are  plotted.  The  value  ■  8.9  corresponds  to  a  charge  density  of  approximately  12.5  x 
-5  3 

10  coul/m  and  a  channel  height  equal  to  d  ■  1*25  Inches.  The  mean  value  of  the  flow 
velocity  is; 


Using  Equation  (98): 


Nr,  I 


tonh  -yUpi,  i"l 


(142) 


(143) 


which  is  the  value  for  ordinary  flow.  The  friction  factor  expression.  Equation  (142),  will 
be  displayed  later  in  comparison  with  the  other  cases. 
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CASE  11.  PARALLEL  PLATES,  VARIABLE 

This  Is  a  case  In  which  the  charge  density  is  assumed  to  vary  in  accordance  with  Equa¬ 
tion  (117).  A  flat  ion  source  is  assumed  to  lie  on  the  y  =  o  plane,  coplanar  with  the  walls, 
as  shown  in  Figure  5.  In  this  idealized  case,  it  is  assumed  that  the  planar  ion  source  does 
nor.  Impede  the  flow  at  the  centerline. 


d*  W 


d  t; 


- •'■fl  P 


W  : 


dp 


Utilizing  Equation  (112); 


for  y  >  y„ 


«  a  [  — — -  1  *  (•»))"* 

'fj  I-  2KA  h  J  '  '  ' 


(144) 


(145) 


Pc„ 


>i~  f  «  >’W  1 

/i.K  I  2h  J 


<  Ju 


v~r 


Let 


Equation  (144)  becomes: 


L  2h 


d*  W  c  w 

- r  -  o  — u 


dp 


41) 

The  solution,  following  Reference  10,  becomes: 


M  '•m 


(146) 


(147) 


(148) 


(149) 


where 


X,  »  o  1)^* 

<9 


(150) 


X,)l8  the  imaginary  Bessel  function.  The  boundary  conditions  are:  (a)  that  the 
solution  must  remain  finite,  and  (b)  at: 


17  *  I,  w  »  0. 

The  integrals  in  this  case  are  not  evaluated  further,  since  this  case  is  not  of  direct 
importance  to  the  evaluation  of  the  test  data.  The  solution.  Equation  (149),  does  serve  to 
Illustrate  the  complexity  of  the  form,  however',  even  when  simplifying  assumptions  are 
used. 
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CASE  III.  ROUND  TUBE,  p  CONSTANT 

'  c 

In  this  case,  round  tube  geometry  Is  considered,  assuming  constant  space  charge,  p^  . 
It  Is  assumed  that  the  wire  does  not  Influence  the  flow;  hence  at  r  *  0,  w  #  0 .  Starting 
with  Equations  (75)  and  (76): 


4  W  1  dW  .1. 

* 

dp* 

(76) 

d rj*  Tj  dr)  °c 

For  constant  charge  density,  p  : 

(151) 

Equation  (76)  becomes: 

4-  -i-  -Np  VI  ^ 

dr)*  r)  dr) 

R« 

dp 

dt 

(152) 

With  boundary  conditions  that  permit  the  solution  to  remain  finite,  and  ar 


r)  *  o  ,  w  *  0. 


The  solution  of  this  equation  is: 


W 


R* 


V)  1 


(153) 


Th's  equation  Is  similar  to  that  for  the  parallel  plate  case  and  the  Hartmann  flow,  except 
that  the  hyperbolic  functions  are  replaced  by  the  Bessel  functions.  Figure  b  Illustrates 
the  flattening  of  the  profile  due  to  the  presence  of  charge  In  the  stream.  N-  •  8.9  repre- 

3 

sente  a  chargo  density  of  approximately  12.5  coul/m  for  the  1.25-lnch-dlameter  tjbe  and 
a  current  of  approximately  580  microamperes  un^r  the  test  conditions. 


The  meat.'  value  of  the  flow  velocity  Is: 


R*  dp  I  r  ^1  ( )  1 


(154) 


As  —  0 ,  Equations  (153)  and  (154)  become  identical  to  Equations  (89)  and  (91)  for  flow 
without  Ionization,  where  p  »  for  zero  Ionization. 


Using  Equation  (97): 


*0  • 


r  '^Pt. 
1  j 


(155) 


where  I,  andi,  are  imaginary  Bessel  functions  (Ref.  11).  Equation  (155)  Indicaioo  that  Ip 
is  S'  complex  function  of  the  charge  density.  It  can  be  shwn  that  as 
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The  variation  of  Iq  will  be  presented  later  In  this  report. 


CASE  IV.  CONCENTRIC  TUBE,  CONSTANT 

In  this  case,  round  tube  geometry  Is  considered,  assuming  constant  space  charge  and 
that  the  wire  does  affect  the  How.  The  boundary  conditions  are: 


w  *  0  ,  Ot  7J  s  I 

and  the  three  different  conditions  near  the  wire  are: 

Case  IVa.  No  charge  present: 

W  s  0 ,  at  7)^  *  0.  004. 

Case  IVb.  Based  upon  Equation  (114): 

W  *  0,  at  »  0.  005. 

Case  IVc.  Based  upon  the  assumption  that  the  corona  region  may  be  of  the  order  of  10 
wire  diameters: 


W  «  0 ,  ot 

The  equation  in  this  case  is  the  same  as 

<t*  W  ^  J_  dW 
d  •?;*  ij  dij 

Case  IVa.  The  solutions  for  this  case 
(103),  and  (104). 


Vo  *  0  0^ 

fo.'  Case  III,  Equation  (152): 
Np  w  •  — —  (-^) 


(152) 


without  space  charge  are  given  in  Equations  (102), 


Case  IVb.  Where  -  0.005: 
'0 


)-«•“*  ]  <*“) 

2R*  I  r  I  _  '  ,  _ _  . 

-—IT  ^  1-^  -  \  ^  *1  ('^  ) 

*  “■  (AT-)  +<>'•*'  (-yffl  )  I,  (0.005  ) 


-  0.204 


/  0.009  \  « 

)  1 


(157) 
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The  friction  factor  can  be  obtained  from: 


h.\ 
I  M  I 


Nr, 


(97) 


where  Is  obtained  from  Equation  (157)  and  and  K,  are  Imaginary  Bessel  functions 
of  the  second  kind. 


Case  IVc.  Where  rj  *  0.04: 


W 


[  I  -  0.1436  I6(,/N^  )  -  0.403  K, 


(158) 


±fr  ,  dp 
h 


(■5r)v 


4-0  403  *  K,  (-/N^)  4-  0.1436  (  )  I,  (o.04/?i^  ) 

c  c  C 


-0.403  (^^=)  K,  (0.04  yfi^)] 

The  friction  factor  can  be  calculated  from; 

,  .  _2_  .  _L 

®  N|H^  \  ft.  dt  '  •«! 

To  obtain  ,  substitute  from  Equation  (157)  into  Equation  (97). 


(159) 


(97) 


Figure  9  Illustrates  the  variation  of  the  profiles  with  a  change  in  for  Cases  IVb  and 
IVc.  Figure  10  illustrates  Case  IVa  and  Case  IVc.  The  velocity  profiles  have  been  normalized 
by  dividing  the  values  by  the  maximum  velocity.  The  effect  of  the  charge  con  be  seen  by  com¬ 
paring  the  two  curves  for  »  0.04  . 


CASE  V.  ROUND  TUBE,  VARIABLE 

In  this  case,  the  round  tube  geometry  is  considered,  and  the  space  charge  is  assumed  to 
vary  with  radius  according  to  Equation  (119)  or  (132).  Equation  (76)  is  applicable  to  this 


case: 

4*  w  1  1 
47)*  V 

dW 

dT) 

a  *•  W  »  Np 

(76) 

From  Equation  (132): 

Pc"  \ 

II 

•V. 
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r,  -- 


4  rg 


Equation  (76)  then  becomes: 


d*  W 


dW 


*  V  itr< 


w 


>'/* 


d  *  V  <ir) 

The  boundary  conditions  for  this  case  are  chosen  as: 

~  O,  at  7)^  I 


(160) 


and  w  remains  finite  at  the  origin. 


A  review  of  the  test  data  on  velocity  profiles  In  Reference  1  indicates  that  when  Ion 
current  flows,  the  profiles  give  the  appearance  that  the  flow  may  not  stick  at  the  central 
wire.  The  assumptions  used  in  this  case,  for  this  reason,  may  be  fairly  realistic.  The 
solution  of  Equation  (160)  is  given  as: 


J.  y'/t 


*  ■  ^  r,  (i-  r!'-)  /•S’'- 1.  ,04 


•.(f  r.'*) 


A  v'* 


-  Ko(C)  Io(€»  ) 

where 

i  .  f  -  y;/«  7,^* 


(161) 

(162) 


It  should  be  noted  that  although  Kq  iCl-'OO  as  C~*'0  •  functions  involved  in  Equation  (161) 
are  well  behaved  and  reach  a  limit  asC  ~*0. 


A  series  of  functions  of  (  Is  presented  In  Appendix  II  to  aid  in  computing  the  velocity 
profiles.  Figure  11  Illustrates  the  effect  of  the  variable  charge  density  on  the  profiles  for 
various  values  of  average  charge  density.  The  peak  value  of  the  curves  decreases  with 
increasing  charge  because  the  curves  axe  computed  for  a  constant  pressure  drop.  As 
the  charge  density  Increases,  the  mean  flow  velocity  decreases.  The  flattening  erf  the 
profile  is  evident.  Figure  12  presents  the  data  of  Figure  11  noimallzed  to  the  value  at 
the  centerline.  In  order  to  determine  tlie  mean  velocity,  the  velocity  profile  curves  were 
integrated.  Cons^uently.  only  specific  values  of  the  average  velocity  were  obtained.  As 
the  charge  density.  Pg  ,  goes  to  zero.  Equation  (161)  becomes  Identical  to  the  form  given 
in  Equation  (89)  or  (91)  for  ordinary  flow  in  a  round  tube  without  charge  effects. 


N 


Pc 


m 


~  0.56 

»  0. 1095  -4^ 

H-  01 


(163) 
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5.55 


>e 


=  0.0598 


_R^ 


(164) 


=  .2.4 

*m  *  -  0.0399  (165) 

Np  =  176 

ooso,-^(-|f-)  (•«) 

The  friction  factors  can  be  calciJated  from  these  equations  for  and  from  Equation  (97) 
for  the  round  tube  friction  factor.  The  variation  of  with  charge  density  will  be  presented 
in  a  graph  in  a  subsequent  section  of  this  report. 


CASE  VI.  CONCENTRIC  TUBE,  /»  VARIABLE 

In  this  case,  the  round  tube  geometry  is  considered,  and  the  space  charge  is  assumed 
to  vary  with  radius  according  to  Equation  (119)  or  (132).  It  is  assumed  that  the  velocity 
goes  to  zero  in  the  vicinity  of  the  wire.  For  purpose  of  analysis,  the  boundary  condition 
is  taken  at  'c  '  'w  •  same  equation  holds  for  this  case  as  for  Case  V. 

d*  W  ^  _l_  dW  W  ^ 

dlj*  7)  d7)  fJ.Wtn 

With  boundary  conditions  of: 

ijj  *  0.04,  W  *0 

■7*1,  w  s  0  . 

the  solution  is  given  by: 

W  .  A,  Io(0  +  B,  K,(e»+  Kolfidf 

-K„(C)  /e*''*  K^lO  d^  ]  (167) 

where  the  coefficients  A  and  B  are  determined  by  application  of  the  two  boundary  conditions: 


A|*  ^1 


B,  *  /9, 


0| 

OrIo(t  ^0^*)-  °0  io(-T  >'•’'*) 


where 


0,  '  —  {-f-) 


(168) 

(169) 

(170) 
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0,  =  lo (-§-  il''*) '<o(t  lo  (-f-  y Ko  (-|-  y/*  )  (i7i) 


T  y^  « 


-Mf  K„(Ode  (172) 


Or  =  *<0  ^-r 

o 


-  ^o{-T-  y,''*)  C*''*  «f)  df  (173) 

0 


The  value  of  is  chosen  as  desired  If  ■  0.04,  as  indicated  In  the  boundary  conditions, 
this  is  inserted  and  the  coefficients  are  determined 


For 

1)  *  0.  04  , 

y,  s  13.17,  *  (7.  55 

W  «  -^(4^)  [-1.(414  I.({)  +  I.  1298  Kflif) 

'  dz.  /  I 

+  I„  (e)/C*/»Ko(€)dC-KoCO/C^»IolC»dC  (174) 

Figure  13  Illustrates  the  velocltjr  profile  associated  with  Case  VI.  The  assumed  value 
for  7]^  has  a  marked  influence  on  the  shape  of  the  profile  near  the  centerline.  It  Is  believed 
that  the  assumption  of  rg  *  10  r*  is  far  too  drastic,  and  that  the  value  ot  Vo"  0.005,  as 
was  used  in  Case  IVb,  may  be  somewhat  more  reasonable.  The  general  influence  of  such 
a  change  In  can  be  obtained  from  the  curves  In  Figure  9  for  the  p,.  ■  constant  case. 
Because  of  the  large  amount  of  work  necessary  for  Case  VI  solutions,  profiles  for 
Vf,  ■  0,005  were  not  calculated. 

The  mean  velocity  for  the  flow  in  Case  IV  with  Vo  "  0.004,  can  be  obtained  by  Inte¬ 
grating  the  velocity  profile.  Thn  friction  factor  can  be  obtained  by  using  this  mean  value 
in  Equation  (97). 

OTHER  CASES 


Two  other  possible  cases  were  studied  in  which  the  charge  density  varied  with  radius. 
The  two  variations  considered  are: 


Pc 


S 


\  RK  / 


(175) 
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The  differential  equation  for  the  density  variation.  Equation  (175),  Is; 


d*  w  (  ^  _3_  R*  ,  dp  » 

d7)*  7)  dT)  4  7)  *  MWm*  ' 


The  differential  equation  for  the  case  of  Equation  (176)  is: 


_w  ^  w  ^  R*  /  dp  V 

7*  V  4  ’?*  '  dz  ' 


An  elementary  solution  is  possible  for  Equation  (177): 


(77) 

[4-  ] 


When  pf.~o.  Equation  (178)  reduces  to  that  for  ordinary  Hagen-Polseuelle  flow.  The 
shape  of  the  profile  Is  given  in  Figure  14. 
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CORRELATION  AND  DISCUSSION 
FRICTION  FACTOR  CURVES 

The  hypothesis  presented  In  the  previous  section  of  this  report  puts  forth  the  concept 
that  the  Ions  In  the  flow  tend  to  set  up  an  Internal  field  within  the  flow.  This  internal  field 
acts  on  the  Ions  to  protide  a  reslstinit  body  force  within  the  flow.  Rased  upon  that  hypothe¬ 
sis,  two  general  forms  for  the  charge-induced  body  force  result.  One  Is  a  mean  value  of 
the  force  given  by: 


F  .  ")m  (59) 

*m  K 

Using  this  average  value,  a  variation  of  friction  factor  for  a  round  pipe  Is  given  by: 

The  second  form  for  the  body  force  Is: 


Using  this  form,  we  found  several  solutions  for  In  the  previous  section. 

The  curves  of  Iq  using  the  mean  value  of  the  body  force  are  presented  in  Figure  15. 
The  equations  for  Iq  for  both  flat  plates  and  round  tubes  Indicate  that  the  friction  factor 
should  Increase  linearly  with  .  Since  Nyo  is  proportional  to  the  square  root  of  the 
corona  current,  then  the  Increase  in  fp  should  occur  as 

Curves  for  the  variation  of  fp  based  upon  the  vari^ible  body  force.  Equation  (179),  are 
shown  In  Figure  16.  This  figure  includes  almost  all  the  cases  considered  in  the  previous 
section.  As  can  be  seen  from  both  Figures  15  and  16,  the  theoretical  predictions  show  the 
increase  In  friction  factor  to  be  much  greater  for  the  flat  plate  case  than  for  the  round 
rube  cases.  The  round  tube  cases  indicate  that  the  effect  of  the  variable  charge  density 
in  a  given  tube  configuration  is  to  increase  the  rate  of  rise  above  the  zero  value.  The 
effect  of  Including  the  influence  of  the  central  wire  or  corona  region  on  the  flow  is  to  dis¬ 
place  the  curves  upward  by  approximately  the  amount  of  the  zero  current  Increase.  Thus, 
the  central  boundary  condition  does  not  have  a  controlling  influence  on  the  pressure  rise 
with  current.  This  same  point  can  be  gleaned  from  rly»  velocity  profile  data  as  well.  Of 
significant  Imponance  Is  that  all  the  theoretical  lines  indicate  an  approximately  linear 
rise  In  the  region  of  under  consideration,  although  the  curves  are  actually  not  linear. 
This  Indicates  that  the  more  complex  equations  also  predict  that  the  Increase  in  Iq  should 
be  prc^rtlonal  to  the  square  root  of  the  corona  current  in  the  region  of  interest. 

The  analysis  of  the  actions  of  space  charge  utilized  the  simplifying  assumption  that 
for  regions  away  from  the  wire: 


Consequently,  was  dropped  in  Equation  (1 15)  for  the  charge  density  in  the  coaxial  case. 
The  analysis  of  Case  V.  which  ignores  the  presence  of  the  wire  but  uses  the  variable  charge 
density  assumption, may,  consequently,  be  somewhat  inaccurate  in  the  region  around  . 
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COMPARISON  OF  VELOCITY  PROFILES 


The  velocity  profiles  In  the  flow  stream  vary  depending  0:1  the  assumptions  made  on 
the  charge  density  and  body  force  term.  In  the  case  of  the  mean  value  force,  the  velocity 
profiles  given  by  Equation  (82)  or  (89)  are  the  same  as  for  ordinary  Hagen-Polseuelle 
flow.  The  profiles  remain  parabolas,  and  the  charge  acts  only  as  a  retarding  Influence  on 
the  flow.  The  action  of  the  charge  In  this  case  can  be  considered  as  effectively  Increasing 
the  viscosity.  This  approach  was  presented  in  Reference  1.  From  Equations  (94)  and  (95), 
the  effective  electroviscosity  could  be  given  by: 


Flat  Plates 


Me  '  M  ('  + 


(180) 


Round  Tubes 

Me  -M 


(181) 


Figure  17  presents  a  summary  of  the  various  velocity  profiles  for  the  round  tube, 
assuming  that  the  body  force  is  variable  in  accordance  with  Equation  (179).  It  can  be 
seen  for  all  the  cases  considered  that,  for  a  fixed  mean  flow  velocity  of  1  ft/sec,  the 
profiles  are  severely  distorted  and  somewhat  flanened  by  the  presence  of  the  ionization. 

A  strong  tendency  exists  for  the  elimination  of  the  large  central  velocitltes  found  In 
ordinary  Hagen-Polseuelle  flow.  A  steepening  of  the  profiles  at  the  walls  of  the  tube  can 
also  be  observed.  This  steepening  should  be  expected  to  result  in  increased  shear  stress 
at  the  wall  with  a  corresponding  increase  in  pressure  drop  in  the  flow.  Such  an  increase 
in  pressure  drop  Is  what  actually  was  observed. 

Figure  18  presents  typical  profiles  for  a  mean  flow  velocity  of  1  ft/sec  in  terms  of  the 
pressure  head  which  corresponds  to  the  local  velocity.  This  presentation  affords  a  better 
comparison  with  the  test  data.  The  test  data  are  given  in  Figures  2  and  19  through  24 
(reproduced  from  Reference  1).  A  comparison  of  the  values  indicates  the  following:  The 
test  data  indicate  that  without  applied  field,  the  presence  of  the  wire  definitely  affects  the 
flow  at  the  lower  Reynolds  numbers.  With  ionization  present,  however,  the  presence  of 
the  wire  does  not  seem  to  affect  the  profiles  as  muclk  In  some  cases  (Fig<jres  19,  23,  and 
24),  the  test  data  indicate  that  the  velocity  across  the  pipe  is  relatively  uniform.  Other 
runs,  plotted  in  Figures  20,  21,  and  22,  indicate  distinct  variations  of  the  profile  in  which 
the  velocity  tends  to  increase  away  from  the  centerline  toward  the  wall,  llie  theoretical 
curves  of  Figure  18  indicate  a  general  flattening  of  the  profiles.  At  higher  ,  turbulence 

could  be  present  and  could  contribute  to  the  flattening  found  in  the  test  data. 

The  Case  V  curves  (Figures  17  and  18)  indicate  a  trend  of  the  profile  to  Increase  from 
the  center  of  the  pipe  outward,  even  though  this  was  not  imposed  as  a  boundary  condition 
on  the  flow.  As  the  ionization  is  increased,  this  gradient  from  the  center  outward  in¬ 
creases.  A  similar  profile  shape  can  be  seen  in  Figures  20  and  21.  The  profiles  calcu¬ 
lated  upon  the  condition  that  the  flow  velocity  near  or  at  the  wire  goes  to  zero  do  not 
seem  to  be  in  as  good  agreement  with  the  test  data  as  are  the  cases  where  the  wire  effect 
is  neglected  in  the  analysis.  No  explanation  for  this  is  available  at  present.  It  would  appear 
acceptable  as  a  first  approximation  to  neglect  the  presence  of  the  fine  wire  and  treat  the 
flow  proUem  as  though  the  wire  did  not  exist.  In  a  general  way,  therefore,  the  test  data 
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on  profiles  at  the  low  Reynolds  numbers  verify  the  theoretical  predictions  that  substantial 
distortions  will  occur  with  the  presence  of  ions. 

COMPARISONS  OF  FRICTION  FACTOR 

Comparison  of  the  theoretical  predictions  for  fo  and  the  test  data  are  presented  in 
Figures  25  through  29.  Figure  25  presents  a  comparison  of  the  value  for  theory  and  for 
air.  These  data  are  based  upon  the  exploratory  test  sequence  described  in  Reference  1. 
The  exploratory  test  sequence  utilized  a  very  short  starting  length,  and  consequently  the 
friction  factors  for  the  case  of  zero  charge  are  higher  than  the  theoretical 


When  this  fact  is  considered,  it  can  be  seen  that  the  data  tend  to  show  a  linear  increase  in 
fQ  with  at  lower  values  of  currents.  Some  falling  off  of  (q  occurs  at  the  higher  values 
of  current,  which  may  possibly  be  due  to  localized  discharge.  The  general  agreement  of 
slope  and  linearity  is  considered  good.  Only  the  curves  for  mean  value  force  and  Case  III 
are  shown  for  clarity. 

Figure  26  presents  a  comparison  of  theoretical  data  with  test  data  for  air  using  the 
longer  starting  length  configuration  described  in  Reference  1,  The  data  for  all  Reynolds 
numbers  fall  on  a  straight  line.  Although  the  magnitude  of  the  increase  is  below  that  pre¬ 
dicted  by  all  the  theories,  the  linearity  achieved  is  encouraging.  The  value 

K  *2.0  X  10'*  m* /vol  f -»«c  , 

often  used  for  air,  is  used  for  both  Figures  25  and  26  in  the  computation  of  charge  number. 

Figure  27  presents  a  comparison  of  theoretical  data  with  test  data  for  oxygen.  Data  for 
both  positive  and  negative  corona  from  an  exploratory  sequence  are  shown.  The  data  indicate 
an  approximately  linear  characteristic  over  a  wider  range  of  charge  number. 

Figure  28  presents  data  for  carbon  dioxide  for  both  positive  and  negative  corona.  Linear¬ 
ity  is  very  good  over  a  very  wide  range  of  charge  number.  Significant  differences  of  slope 
appear  between  the  negative  and  positive  conditions. 

Figure  29  illustrates  the  effect  with  nitrogen.  These  data  are  very  limited  and  are  only 
indicative  of  the  trend.  Lines  were  drawn  only  to  relate  the  points,  and  are  not  meant  to 
indicate  that  a  linear  relationship  necessarily  exists. 

The  ccxnparisons  indicated  in  Figures  25  through  29  are  made  with  the  curves  for  Case 
III  and  the  mean  value  curve  for  the  round  tube.  For  a  comparison  with  the  other  theoretical 
cases.  Figure  16  should  be  consulted.  In  most  of  the  cases  presented,  the  data  indicate  an 
approximate  linear  characteristic,  at  least  in  the  lower  ranges  of  corona  current.  Test  data, 
however,  generally  fall  below  the  theoretical  lines. 

EFFECTS  OF  MOBILITY  AND  NONUNIFORMITY  OF  CHARGE  DENSITY 

A  possible  explanation  for  the  observation  that  test  data  fall  below  the  theoretical  values 
may  lie  in  the  value  of  mobility  and  in  the  degree  of  uniformity  of  ionization  in  the  tube. 

The  value  of  mobility  affects  the  charge  number  as: 
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\  “  • 

A  dtnall  change  in  the  value  of  mobility  can  lead  to  substantial  changes  in  the  magnitude 
of  the  charge  number.  Since  precise  values  of  mobility  are  difficult  to  obtain  and  since 
the  values  vary  depending  upon  the  contaminants  and  impurities  in  the  gas  stream,  it 
would  not  be  surprising  to  find  discrepancies  in  the  charge  density.  As  was  indicated  in 
Reference  1,  the  moisture  content  of  the  gases  was  relatively  high.  In  the  case  of  air, 
the  moisture  content  was  approximately  1%,  even  though  the  gas  was  commercially  dried 
air.  The  presence  of  water  vapor  could  lead  to  ions  of  lowered  mobility.  If  the  average 
mobility  was  indeed  lower,  then  the  charge  densities  calculated  using  the  K  values  given 
in  Table  1  could  be  somewhat  inaccurate. 

In  the  calculation  of  charge  density,  it  was  tacitly  assumed  that  during  the  corona  dis¬ 
charge  the  ionization  would  be  uniform  along  the  wire  and  remain  uniform  as  the  applied 
field  and  corona  current  increased.  This  is  an  idealized  assumption  which  probably  can¬ 
not  be  achieved  in  practice.  It  is  the  characteristic  of  corona  discharge  to  tend  to  local¬ 
ize  in  certain  discrete  points  as  the  field  increases.  During  corona  discharge  in  an  open 
room,  the  characteristic  hissing  emenating  from  these  local  points  can  be  heard  quite 
readily.  During  such  testing  in  the  dark,  such  local  points  were  observed  at  several  lo¬ 
cations  along  the  electrode  as  pinpoints  of  light  surrounded  by  a  violet  glow.  At  these 
points,  the  current  may  be  considerably  larger  than  the  average  for  the  wire.  Such  local 
dischargee  would  result  in  local  high  charge  densities  and  a  nonuniform  distribution  of 
ions  along  the  tube.  Thus,  a  substantial  portion  of  the  total  charge  might  be  located  in 
discrete  narrow  paths  that  lay  radially  between  the  electrodes.  Under  this  condition,  cal¬ 
culating  the  average  charge  density  from  the  measured  corona  current  could  lead  to  a 
substantial  overestimation.  Such  possible  nonuniformities  in  charge  density  could  have 
resulted  in  values  of  for  Individual  test  points  that  were  higher  than  actually  existed 

over  the  major  portion  of  the  tube. 

A  second  action  in  the  test  could  lead  to  nonuniform  charge  distribution.  There  is  a 
possibility  of  slight  deformation  and  sagging  of  the  wire.  If  the  wire  were  to  move  out  of 
concentricity  so  that  it  is  closer  to  one  side  of  the  channel  than  the  other,  a  more  Intense 
field  would  result  in  the  region  of  close  spacing.  A  greater  local  ionization  could  result 
which  would  tend  to  traverse  the  shortest  lines  of  force  to  the  nearest  wall.  Such  an  action 
could  lead  to  undesirable  nonuniformities.  Prior  to  setting  up  the  tests,  we  made  a  series 
of  calculations  of  the  droop  due  to  gravity  and  e!'>ctrostatlc  forces.  The  wire  tension  was 
adjusted  to  keep  the  wire  straight  within  1/64  Inch  under  the  action  of  gravity.  Strength 
limitations  of  the  wire  did  not  permit  reducing  the  deflections  much  below  this  value.  If 
the  wire  were  not  in  perfect  center,  then  the  electrostatic  forces  could  act  to  move  the 
wire  further  out  of  true  concentricity.  Consequently,  some  nonuniformity  of  space  charge 
may  have  resulted  from  the  deformations  that  can  occur  during  test.  It  is  considered 
doubtful,  however,  that  such  localized  discharges  could  have  contributed  much  to  the  in¬ 
crease  in  pressure  drop  and  profile  distortion.  Consequently,  it  is  believed  that  their 
primary  effect  will  result  in  the  overestimation  of  as  discussed  above.  Such  non- 

uniformities  in  charge  density  may  possibly  explain  why  the  test  data  did  not  show  as 
large  a  rate  of  increase  in  (q  as  was  predicted  by  the  theory. 
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SUMMARY 


In  order  to  seek  an  explanation  of  the  unusually  high  pressure  drops  reported  In  Refer¬ 
ence  1,  we  conducted  an  analysis  based  upon  the  hypothesis  that  the  Ions  In  the  stream 
tended  to  set  up  an  electiical  field  Internal  to  the  stream.  Although  it  Is  not  possible  to 
derive  a  unique  solution  for  the  internal  field  with  precise  mathematics,  it  is  believed  that 
an  expression  of  the  form 


= '  (x) 

is  a  valid  representation  of  the  field.  Based  upon  the  mean  value  approach  for  the  charge 
density,  field,  and  velocity  in  the  channel,  explicit  solutions  were  olMained  which  predicted 
large  increases  in  overall  pressure  drop  when  ionization  is  present  and  is  confined  In  the 
tube  by  the  applied  field.  The  increases  predicted  were  of  the  correct  order  of  magnitude 
as  compa’'ed  to  the  data  of  Reference  1.  The  mean  value  approach  is  considered  to  be 
mathematically  sound.  The  expression 


satisfied  the  electric  field  equations  Identically,  and  is  a  permissible  but  not  a  unique 
representation  for  E*  . 

If  some  ion  slippage  could  occur,  such  that  the  ions  would  not  move  completely  against 
the  stream  velocity,  an  expression  of  the  form 

) 

could  be  utilized,  where  Oq  is  a  constant  which  indicates  the  amount  of  relative  slip  present. 

If  this  constant  is  included  in  the  equations,  then  the  charge  number  is  merely  multiplied 
by  the  constant  and  we  get  Op  .  Thus,  the  effect  of  Including  such  an  assumed  slip  is 

to  reduce  the  effectiveness  of  the  charge  number.  Inclusion  of  such  a  term  will  not  change 

any  of  the  solutions  presented  but  would  change  only  the  magnitude  of  the  effective  charge 

number.  It  would  be  possible  to  obtain  empirical  values  for  Og  by  determining  the  change 

in  effective  charge  number  necessary  to  make  the  slopes  for  the  theoretical  curves  in 

Figures  26,  27,  and  28  agree  with  the  test  data.  No  analytical  representation  is  available 

at  this  time,  however,  to  fully  Justify  the  Inclusion  of  such  a  slippage  constant.  The  fact  ' 

that  theory  tends  to  show  greater  pressure  drops  than  experiment  may  be  an  Indication 

that  some  slippage  actually  does  exist  in  the  stream.  * 

Use  of  the  expression 


in  the  fluid  equations  led  to  a  linear  equation  with  variable  coefficients.  This  equation 
was  solved  for  several  channel  geometries  and  various  representative  charge  distribu¬ 
tions.  The  solutions  for  these  cases  indicated  increases  in  pressure  drop  and  friction 
factor  similar  to  those  found  with  the  mean  value  approach.  The  rate  of  Increase  with 
charge  density  for  these  cases  was  somewhat  higher  than  predicted  by  the  mean  value 
approach.  All  the  cases  indicated  thatlarge  changes  to  the  velocity  profile  would  occur 
in  a  channel  with  bounded  charge  present.  The  correladon  at  the  velocity  profile  test  data 
and  the  tlieoretical  curves  is  considered  to  be  good.  In  some  cases,  fairly  close  represen¬ 
tation  was  achieved. 
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Putting  the  equations  In  nondlmenslonal  form  resulted  In  the  appearance  of  two  new 

nonditnenslonal  parameters.  These  two  parameters  are  the  charge  number,  Nn  .  and 

»  c 

the  charge  Reynolds  number  Np  .  represents  the  ratio  of  the  charge  forces  to  the 

viscous  forces,  and  Np^  represents  the  ratio  of  the  Interla  forces  to  the  charge  forces. 

Appendix  I  presents  a  detailed  analysis  of  these  parameters.  Based  upon  the  kinetic 
tteory.  It  was  shown  that  the  charge  number  is  a  direct  function  of  the  number  charge 
density  ratio  and  the  Inverse  function  of  the  Knudsen  number  squared.  If  any  appreciable 
degree  of  ionization  is  present,  the  Knudsen  number  acts  as  an  extremely  large  amplifi¬ 
cation  factor.  This  term  provides  a  fundamental  explanation  of  the  surprisingly  large 
effects  on  pressure  drop  and  profile  observed. 

The  overall  correlation  achieved  between  the  theoretical  predictions  and  the  test  data 
Is  believed  to  be  gooc'.  The  linearity  achieved  with  vs.  is  considered  to  agree 

with  the  predictions.  A  clear  picture  of  the  linearity  can  also  be  seen  from  the  additional 
data  presented  in  Appendix  III.  The  analysis,  however,  predicted  a  greater  increase  in 
friction  factor  for  all  cases  considered  than  was  actually  observed.  The  factors  which 
can  affect  this  Include  the  mobility  of  ions,  the  uniformity  of  the  charge,  and  the  assump¬ 
tion  that: 


The  mobility  of  ions  is  difficult  to  determine  exactly  for  a  given  situation  and  consequently 
can  give  rise  to  error.  Nonunlformlties  of  corona  discharge  In  the  stream  can  lead  to  an 
overestimation  of  the  charge  density.  The  form  of  the  Induced  field, 

---T 

could  be  an  overestlmatlon  of  the  Influence  of  the  charge.  The  term  j^/p  k  from  Equation 
(52a) 


could  act  to  reduce  the  net  induced  field  in  tlie  stream. 

Following  completion  of  the  work  reported  herein,  it  was  learned  that  a  paper  by 
Stuetzer  on  ‘‘Appare;it  Viscosity  of  a  Charged  Fluid”  might  be  applicable  to  this  prob¬ 
lem.  A  review  of  that  paper  Is  presented  in  Appendix  IV.  Stuetzer  utilized  a  similar 
analytical  approach  that  produced  average  values  Identical  to  the  mean  values  presented 
in  this  report.  This  agreement  is  considered  to  provide  strong  suppon  for  the  concepts 
set  forth  in  this  analysis. 
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CONCLUSIONS 


From  the  analysis  carried  out  In  this  report  and  from  a  comparison  of  the  theoory  with 
the  test  data  of  ASD-TDR-63-842,  the  following  conclusions  can  be  drawn. 

1.  The  analysis  using  the  mean  value  the  body  force  provides  a  solution  which  Is 
sound  from  a  mathematical  standpoint.  The  analysis  Indicated  that  large  Increases  in 
pressure  drop  and  friction  factor  could  be  expected  if  ionization  Is  present  In  the  stream 
and  the  charge  Is  bound  to  the  physical  system.  The  Increases  predicted  by  the  mean 
value  approach  are  of  the  same  order  as  those  actually  observed. 

2.  The  analysis  utilizing  the  first  order  approximation,  that 


indicated  that  large  Increases  In  pressure  drop  and  large  distortions  to  the  velocity  profile 
would  occur  If  a  bound  charge  existed  In  the  stream.  Good  correlation  with  the  test  data  was 
achieved. 

3.  The  predictions  of  theory  that  the  increases  in  friction  factor  and  pressure  should 
vary  directly  as  the  charge  number,  Nn  ,  or  as  the  square  root  of  corona  current  were 

substantiated  by  comparisons  with  the  test  data. 

4.  Two  new  nondimenslonal  parameters  appeared,  the  charge  number,  and  the  charge 
Reynolds  number.  These  two  parameters  control  the  effects  of  charge  on  flow.  The  charge 
number  is  a  function  of  the  square  of  the  Knudsen  number.  The  Knudsen  number  acts  as 
an  amplification  factor  on  the  charge  present,  and  provides  the.  <.ey  to  the  surprisingly 
large  actions  of  charge  on  flow  of  dense  gases. 

5.  The  hypothesis  that  the  Ions  set  up  an  internal  field  acting  to  retard  the  flow  is  con¬ 
sidered  to  be  substantiated  by  the  comparison  of  theoretical  values  with  test  data.  The 
solutions  for  velocity  profile,  however,  may  not  necessarily  be  unique,  and  further  work 
is  needed  in  this  direction.  Stuetzer’s  analysis  on  apparent  viscosity  provided  further 
substantiation  for  the  validity  of  the  Induced  Intemd  field. 

6.  The  net  Impact  of  these  charge  phenomena  depends  upon  further  careful  evaluation 

of  the  charge  effects  In  both  internal  and  external  flow.  It  Is  necessary  to  define  the  limits 
of  the  action  as  dimensions  vary  and  as  flow  velocities  Increase.  Possible  critical  or 
limiting  values  for  and  Np^  must  be  determined.  If  further  research  Indicates  that 

the  charge  phenomenon  is  active  in  laminar  external  flows,  then  we  can  envision  many 
ultimate  practical  applications.  The  phenomenon  could  conceivably  be  applied  to  the 
control  of  laminar  boundary  layers  and  heat  transfer  under  laminar  conditions.  It  would 
appear  useful  as  a  basic  flow  investigative  technique  for  use  In  the  study  m  boundary 
layers.  It  could  possibly  be  used  as  a  flow  trigger  In  regions  of  transition  or  separation 
In  the  gas  stream.  It  is'not  possible  at  this  time  to  delineate  explicitly  the  techniques  of 
Incotporating  the  mechanism  in  actual  operating  systems;  it  Is  believed,  however,  that 
applications  will  evolve  as  experience  is  gained  In  working  with  and  controlling  the  phe¬ 
nomena. 
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Figure  4.  Round  Tube  Channel  with  Central  Wire 


Figure  5.  Parallel  Plate  Flow  Channel 
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Figure  8.  Case  m.  Theoretical  Round  Tube  Velocity  Dlstribidion,  p  :  Constant 
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heoretical  Concentric  Tube  Velocity  Distribution  Case  IVb  and  IVc,  p 
Tube  * 


K  Constant 


Theoretical  Concentric  Tube  Velocity  Distribution  Case  IV,  p^  e  Constant  in  Tube 
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Figure  11.  Case  V.  Variation  of  Velocity  Distribution  at  Constant  Pressure  Drop  and  with 
ot  77  -  '/I 
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Fli(ure  15.  Eti'ec 


54 


tHIIIMtg.l 


EUSH 


Figure  17.  Typical  Velocity  Diatributions  at  w  =1  ft/aec. 
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Figure  18.  Typical  Profiles  in  Terms  of  Pressure 
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Figure  19.  Vulocity  Profile  Without  lonlzotlon.  Concentric  0.00  S  -Inch  Wire,  w  s  2.  35  ft/ sec. 
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Figure  2D.  Velocity  Profile  Data,  N„  >12)0 
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Figure  23.  Velocity  Profile  DalR,  N„  =  1745 
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Ficure  25.  E;slcritory  Teat  Data  Showing  Variation  of  (_  with  Charge  Number  for  Air 
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Figure  26.  Variation  o(  with  Charge  Number  for  Air 
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with  Charge  Number  for  Oxygen 
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Figure  29.  Effect  of  Charge  Number  on  f^^  for  Nltrogeo 
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^  MX 


A>  "mP 


) 


.  J_ 

^  A>  *mK 


(183) 


Tills  represents  the  ratio  of  the  charge  forces  to  the  inertia  forces  In  the  stream.  The 
Inverse  of  Equation  (183)  is: 


(184) 


The  factor  In  the  brackets,  then,  becomes  an  effective  charge  Reynolds  number. 


(185) 


This  interesting  parameter  indicates  that,  as  the  charge  density  disappears,  the  charge 
Reynolds  number  becomes  very  large  and  the  flow  field  Is  dominated  by  the  inertia  forces. 
If  there  is  appreciable  ionization  in  the  stream  and  the  density  of  the  gas  is  low,  then 
charge  effects  could  appear.  If  the  velocity  is  very  low,  charge  effects  could  also  appear. 

If  dimensions  are  very  large  and  velocity  and  density  are  relatively  low,  the  charge  could 
also  influence  the  flow.  The  data  could  be  plotted  against  the  charge  Reynolds  number,  but 
the  plots  of  Iq  vs.  at  constant  corona  current  (constant  charge  density)  and  fo  vs. 

N„  are  believed  to  present  the  data  adequately. 

It  is  of  considerable  interest  to  study  the  flow  charge  number  further: 


Np  *  - ^ 

Pc 


Since  both  ^  and  K  depend  iqxm  the  molecular  actions  in  a  gas,  this  ratio  may  be  expressed 
in  completely  different  terms. 

From  the  kinet.c  theory  of  gases  the  viscosity  of  a  gas  can  be  given  by  (Ref.  13): 

(186) 

where  N  is  the  number  of  molecules  per  tudt  volume, 
m  Is  the  mass  of  a  molecule, 

T  is  the  mean  speed  of  a  molecule,  and 
X  Is  the  mean  free  path. 


The  constant.  1/3,  is  fro.m  the  elementary  derivation  of  viscosity.  More  complete  theories 
modify  this  constant.  To  account  for  this,  let  the  vtscoslty  of  the  gas  be  represented  by: 


70 


RTD-TDR-63-4009 


^  =  a,  N  m  c  X 


(187) 


where  a,  is  a  suitable  constant. 


From  the  kinetic  theorj',  the  value  of  mobility  in  terms  of  molecular  action  can  be  obtained 
(Ref.  2)  by: 


q  X| 


m.  c 


(188) 


where  q  is  the  charge  per  ion. 

Xj  is  the  mean  free  path  of  an  ion, 
iDj  is  the  mass  of  the  ion. 

T  is  the  mean  speed  of  the  neutral  molecules  in  the  gas  through  which  the  ion  is 
moving,  and 

b  is  a  constant  which  depends  upon  the  theory  used. 

Coblne  (Ref.  2)  stKn/s  that  for  Ions  of  the  parent  gas: 

X,  *  2''*X 


(189) 


An  experimental  observation  Is  that  the  electrical  forces  around  an  ion  cause  its  mean 
free  path  to  be  much  shorter  than  that  of  the  neutral  molecule  (Ref.  4).  The  form  for  K 
given  in  Equation  (188)  Is  sufficiently  accurate  for  our  use,  since  the  value  for  the  constant, 
b  ,  can  be  chosen  to  account  for  the  theory  used  and  rhe  mean  free  path.  F  rom  Equations 
(187)  and  (188): 


uKzo.NmcX  Xbq  — ^ 

me 


Using  Equation  (189): 

/*  K  »  z'/*  0,  b  q  N  X* 
but 

Letting  Ob  s  Cj  and  substituting  Equation  (192)  gives: 


Substituting  into  the  charge  number: 


'  t  (^)  (^) 


(190) 

(191) 

(192) 

(193) 

(194) 
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The  Knudsen  number  is  defined  as: 


(195) 


It  is  the  ratio  of  the  mean  free  path  to  a  reprerentatlve  dimension  of  the  flow  system.  Sub¬ 
stituting  Equation  (195)  Into  the  charge  number; 


(196) 


N|  /N  represents  the  number  of  charges  per  neutral  molecule  and  Is  dimensionless.  Cj  , 
a  constant.  Is  dimensionless,  as  is  the  Knudsen  number. 


(197) 


This  equation  Indicates  why  the  lonlratlon  has  such  a  surprisingly  large  effect  on  the  flow. 
If  ionization  Is  present  In  any  appreciable  quantity,  then  the  Influence  of  the  Knudsen  num¬ 
ber  can  be  very  great.  Since  In  ordinary  gas  flows  at  standard  conditions: 


and  using  L  *  .  (i )  cm  : 


X  *  0  1 lO’*  )  cm 


*  s  0  (10*'®  ) 


Considering  the  constant,  cj  ,  based  upon  the  simpler  theory: 

I 


b  *  I. 


Therefore: 


=  o 


#> 

yr 

C|  K 

3 

3 

X  lo' 

/r 

.^*0 

_Nl\ 

N  ! 

10 

N  / 


This  indicates  that  the  charge  number  will  be  of  the  order  of  10^  times  the  charge  ratio. 
This  large  amplification  effect  of  the  Knudsen  number  .ndlcates  why  even  small  charge 
densities  can  have  significant  effects  on  flow.  Since  typiCi  I  charge  densities  In  the  tests 
-10  -9 

were  from  10  to  10  ,  the  magnitude  of  Np  Is  of  the  order  of  1  to  10.  This  value  Is 

“c 

In  line  with  the  values  of  Np  presented  In  the  body  of  the  report. 

Other  imponant  ratios  such  as  those  relating  to  compressibility  (Mach  Number)  or 
licat  transfer  can  also  be  obtained  as  necessary  for  specific  problems. 
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APPENDIX  Ji 


CURVES  FOR  CALCULATING  CASE  V  AND  CASE  VI  VELOCITY  PROFILES 


Universal  curves  for  calculating  the  velocity  profile  for  round  tube  geometry  with  vari¬ 
able  space  charge  (Cases  V  and  VI)  are  presented  in  this  appendix.  For  these  calculations, 
C  Is  a  useful  coordinate  defined  by  Equation  (162): 

f  ■  -i-  y  '4  ^*4 


By  using  these  parameters.  Equation  (161)  for  Case  V,  where  the  wire  does  not  alfect  the 
flow,  and  Equation  (167)  for  Case  VI,  where  the  flow  is  assumed  at  rest  at  or  near  the  wi’-e, 
we  can  use  the  curves  to  compute  the  profiles  for  specific  values  of  corona  current. 
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Figure  II-l  thru  n-lO.  Universal  Curves  for  Calculating  Velocity  Profilea,  Cases  V  and  VI 


Figure  n-2. 
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appendix  III 


PRESSURE  DROP  DATA  AND  MODIFIED  CORONA  WIRE  TEST 
A.  VARIATION  OF  PRESSURE  DROP  WITH  CURRENT 

In  the  body  of  this  report  the  friction  factor,  f ^ ,  was  compared  with  the  charge  number, 

.  Agreement  was  shown  between  the  theory  and  the  test  data.  A  clearer  relationship 

between  the  ion  influences  on  flow  and  the  charge  density  can  be  obtained  by  considering 
the  variation  of  pressure  drqj  with  increasing  corona  current.  Figures  lll-l,  -2,  and  -3, 
present  typical  data  talcen  from  ASD-TDR-63-842.  In  these  figures  the  pressure  drop  is 
plotted  against  the  square  root  of  the  corona  current. 

Figure  III-l  is  based  upon  data  taken  for  the  exploratory  sequence  for  air  with  the  posi¬ 
tive  corona.  At  the  lower  Reynolds  numbers  the  data  points  fall  on  straight  lines  indicating 
good  agreement  with  the  theory.  In  Figure  III-2  data  is  plotted  from  anotlier  test  sequence, 
which  indicates  an  Interesting  phenomenon.  The  linear  correlation  is  good,  but  it  can  be 
noted  that,  at  the  higher  Reynolds  numbers,  a  region  of  very  low  elope  occurs  at  low  current 
values.  As  current  values  increase,  there  is  a  linear  Increase  at  a  greater  slope.  We  inter¬ 
pret  this  phenomenon  to  indicate  that  a  much  larger  charge  density  is  required  at  the  edge 
of  the  transition  zone  to  make  its  presence  felt  in  the  flow  and  that  there  may  be  a  critical 
value  for  the  charge  density  before  the  action  becomes  effective.  In  Figure  III-l,  the  gen¬ 
eral  falling  off  of  the  rate  of  increase  at  the  higher  values  of  current  may  be  due  to  local¬ 
ized  discharge  effects.  This  lower  slope  at  the  high  currents  may  be  due  to  corona  heating 
of  the  stream,  which  would  cause  only  a  slight  Increase  in  pressure  drop.  The  magnitude 
of  the  heating  effects  are  considered  in  Reference  1. 

Figure  III-3  presents  data  for  CO2  with  a  positive  corona.  At  the  lower  flow  rates  the 

linear  characteristic  Is  evident;  at  the  higher  flow  rates,  however,  the  effect  of  charge 
shows  a  pronounced  nonlinear  characteristic.  Figure  III-4  presents  additional  data  for  CO2, 
at  two  widely  separated  Reynolds  numbers  and  a  wider  range  of  current.  At  low  currents 
and  lowNf^  ,  the  data  of  Figure  I1I-4  compare  well  with  Figure  III-3.  At  high  Nj,  ,  the 
ion  effects*are  still  present  but  the  curve  is  nonlinear  throughout  Its  range. 


All  of  the  figures  indicate  excellent  correlation  in  the  laminar  flow  regime  at  the  lower 
values  of  current.  In  this  region  the  theoretical  predictions  are  adequately  substantiated 
by  the  test  data. 

B.  TESTS  WITH  MOEMFIED  CORONA  WIRE  END  CONDITION 

A  brief  retest  of  the  corona  wire  was  made  because  of  considerable  discussion  over  the 
drift  motion  of  the  ions.  It  had  been  suggested  that  the  ions  would  merely  drift  downstream 
slightly  due  to  the  stream  velocity  and  that  a  local  Increase  in  current  at  the  outlet  of  the 
pipe  at  the  outer  electrode  would  occur.  To  determir.e  whether  the  end  conditions  of  the 
corona  wire  and  pipe  were  significant,  an  insulating  tubing  was  inserted  over  the  corona 
wire  and  extended  1-1/4  inches  In  from  the  outlet  of  the  pipe.  Data  were  taken  for  CO2  at 

a  range  of  currents  and  Reynolds  numbers.  No  changes  to  the  charge  phenomenon  reported 
in  ASD-TDR-63-842  were  ^served.  It  Is  concluded  from  this  test  that  the  exact  end  condi¬ 
tions  for  the  corona  discharge  are  not  necessarily  critical  for  the  phenomenon  to  occur. 
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APPENDIX  IV 


COMPARISON  WITH  STUETZER’S  PAPER 

F  ollowing  completion  of  the  body  of  this  report,  the  author  became  aware  of  the  possible 
similarity  Lietv.'een  the  plienomenon  dljcussed  by  Stuetzer  in  his  paper,  “Apparent  Viscosity 
of  a  Charged  Fluid  "  Reference  12,  and  the  phenomenon  discussed  herein.  Stuetzer,  in  his 
work  with  ionized  gases,  discusses  the  apparent  viscosity  of  a  charged  fluid.  The  Initial 
stages  of  .Stuetzer's  analysis  are  remarkably  similar  to  those  of  the  analysis  presented  in 
this  repon,  although  the  two  research  efforts  were  conceived  and  conducted  completely 
independent  of  each  other. 

Stuetzer  starts  with  the  same  Iwsic  electric  field  and  fluid  equations*  and  places  a  condi¬ 
tion  on  the  current  that: 


‘2  '  /»c  ‘ 

-h 

The  current,  in  Equation  (198),  is  the  net  current  in  the  axial  direction.  If  Ij  is  not  zero, 
there  is  a  net  charge  flow  down  the  tube,  and  a  corresponding  “corona-wind"  pressure  rise 
can  occur  along  the  pipe  due  to  this  charge  motion. 


Using  Equation  (49): 


j'  s  (u^  d-  KE*  ) 


and  the  mean  velocity  defined  by: 


(49) 


I  r  . 
*  TT  /  " 


-h 


and  substituting  into  Equation  (11): 


.  dp 

u  7  •  -  -r —  +  p  E,  =  0 

^  dt  ^ 


Stuetzer  obtains; 


dy* 


^  d  w  +  * 


!  »  I 

i-c  I 


dx  ~  Z^jIk  ZhK 


j  w  dy 


«  0 


— h 


(199) 


(H) 


(200) 


where  |p^|  is  the  absolute  value  of  the  charge  density.  Integrated,  this  becomes  Stuetzer’s 
Equation  (11): 


(201) 


h*/-^  - 

n  1 

\  dx 

2h  /  K  • 

V  m  MM  A 

\ 

5  +  />c 

3  K  / 

The  symbols  used  in  Stuetzer’s  paper  have  been  changed  to  conform  to  those  used  here. 
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The  term  containing  is  the  ‘‘corona  wind”  pressure  rise  term.  For  =  0 ,  this  expres¬ 
sion  becomes  identical  to  Equation  (94)  for  parallel  plates. 

Likewise,  Stuetzer  derives  two  expressions  for  apparent  viscosity  for  flat  plates  and 
round  tubes: 


Flat  Plates 


Round  Tubes 


{\Pc\ 


^E  '  M  +(  l^cl  -^  ) 

These  equations  are  Identical  to  the  mean  value  approximations  given  by  Equations  (180) 
and  (181).  Stuetzer  also  calls  attention  to  the  similarity  of  this  type  of  flow  to  Hartmann 
flow  of  MHD  and  draws  an  Interesting  comparison  with  the  effective  viscosity  Increase 
in  both  types  of  flows. 

The  analysis  given  in  the  body  of  this  report  using  the  mean  value  of  the  body  force 
leads  to  the  same  value  of  apparent  viscosity  as  Stuetzer  obtained  by  his  method.  The  two 
approaches  are  similar,  but  the  tedrjques  used  to  Introduce  the  average  value  are  different. 
In  discussing  limitations  on  the  effect  of  the  Ions  on  the  flow,  Stuetzer  Indicates  that  the 
time  required  for  the  fluid  to  flow  down  the  length  of  the  pipe  must  be  short  relative  to  the 
time  required  for  the  charge  to  change  Its  lateral  distribution.  In  the  initial  hypothesis, 
for  typical  values  of  mobility  and  field.  Equation  (1)  Indicated  a  cross  channel  mean  speed 
of  the  order  of  100  m/sec,  while  the  speed  down  the  pipe  was  of  the  order  of  1  m/sec. 
Stuetzer’s  criteria  are: 

Transit  time  down  the  channel,  , 


'■  -(f  T^) 


Transit  time  across  the  channel,  ty  , 


ty  »  O 


(-fir) 


Stuetzer  states  that  his  solution  is  valid  only  when  ty  » ty  .  Examining  the  transit  time 
across  the  channel  further,  consider  Equations  (2),  (3),  and  (30). 
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Therefore: 

—  =  lA 

dy  '  < 


Assuming  a  constant  space  charge: 


Let: 


Let: 


V 


( 


4-  C 


w  =  O,  ot  y  =  0 


£y 

y 


dt 


y  +  c  . 


(207) 

(208) 

(209) 

(210) 
(211) 
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For  the  case  discussed  above: 

»i  * 


Since  ty  «  .  it  appears  that  Stuetzer's  criteria  may  not  be  fully  applicable.  One  possible 

explanation  is  that  in  the  region  between  the  sheaths  occurring  at  each  electrode,  the  local 
field  may  be  low  enough  that  the  local  transverse  ion  drift  velocity  is  very  low  and  the 
transit  time  in  this  region  may  be  significantly  longer. 

From  an  overall  comparison  of  Stuetzer's  paper  and  this  report,  it  appears  that  the  mean 
value  approach  is  similar  to  that  presented  by  Stuetzer.  The  broader  interpretation  used  in 
this  study,  that: 


-0  ,  ol  y  *  0 


0  1 1  s«c  ) 

0  (  10**  t*.’  ) 


(212) 


or 
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was  not  considered  by  Stuetzer.  Corsequently,  he  limits  his  work  to  situations  where  the 
velocity  profiles  do  not  change,  whereas  our  test  results  indicated  severe  distortions  of 
the  velocity  profiles.  The  Stuetzer  paper  is  an  excellent  work,  however,  and  the  author 
regrets  not  recognizing  its  applicability  to  the  charge  phenomenon  sooner  since  it  adds 
solid  substantiation  to  the  hypothesis  used  in  this  report. 
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